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Abstract 

It has been empirically observed in both experiments and particle-in-cell simula- 
tions that space-charge-dominated beams suffer strong growth in statistical phase- 
space area (degraded quality) and particle losses in alternating gradient quadrupole 
transport channels when the undepressed phase advance ao increases beyond about 
85° per lattice period. Although this criterion has been used extensively in practical 
designs of strong focusing intense beam transport lattices, the origin of the limit has 
not been understood. We propose a mechanism for the transport limit resulting from 
classes of halo particle resonances near the core of the beam that allow near-edge 
particles to rapidly increase in oscillation amplitude when the space-charge intensity 
and the flutter of the matched beam envelope are both sufficiently large. When cou- 
pled with a diffuse beam edge and/or perturbations internal to the beam core that 
can drive particles outside the edge, this mechanism can result in large and rapid 
halo-driven increases in the statistical phase-space area of the beam, lost particles, 
and degraded transport. A core-particle model is applied to parametrically analyze 
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this process. Extensive self-consistent particle in cell simulations are employed to 
better quantify space-charge limit and verify core-particle model predictions. 
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1 Introduction 

The maximum transportable current density of an ion beam with high space- 
charge intensity propagating in a periodic focusing lattice is a problem of 
practical importance[l,2]. Accelerator applications such as Heavy Ion Fu- 
sion (HIF), High Energy Density Physics (HEDP), and transmutation of 
nuclear waste demand a large flux of particles on target. A limit to the 
maximum current density can result from a variety of factors: instability of 
low-order moments of the beam describing the centroid and envelope, insta- 
bility of higher order collective modes internal to the beam, growth in sta- 
tistical phase-space area (rms emittance growth), excessive halo generation, 
and species contamination associated with issues such as the electron cloud 
problem. Simulations were first used to analyze the maximum current density 
transportable in quadrupole channels[3,4] and provided guidance beyond ini- 
tial heuristic estimates [5]. Experiments later obtained results consistent with 
simulations [1,2]. 

The present work describes a promising new approach toward predicting the 
maximum transportable current density in a periodic quadrupole lattice due 
to intrinsic space-charge limits [6]. Previous studies to predict space-charge 
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related transport limits in the absence of focusing errors and species contam- 
ination have not proved fully successful beyond a moment level description of 
low-order beam instabilities. Although moment-based centroid and envelope 
descriptions reliably predict regions of parametric instability where machines 
cannot operate[8,7], such models are overly optimistic when compared to sim- 
ulations and experiments which observe degraded transport due to emittance 
growth and particle losses where the moment models predict stability[l,2,3,4]. 
On the other hand, higher-order collective mode theories based on the equilib- 
rium KV distribution[10] predict broad parametric regions of instability where 
stability is observed in simulations with more realistic distributions [3, 4] and 
in experiment [1,2]. The space-charge limit model proposed is based on parti- 
cles oscillating outside, but near the beam edge exchanging energy with the 
oscillating space-charge field of a envelope matched beam core leading to in- 
creased particle oscillation amplitude, emittance blow up, and particle losses. 
This model can be applied to a wide range of matched core distributions and 
does not require an equilibrium core - which circumvents the practical problem 
of no smooth core equilibrium distribution being known. The increased un- 
derstanding the origin of the observed limits obtained promises more reliable 
design of optimal intense beam transport channels. 

We denote the phase advance of particles oscillating in a periodic focusing 
lattice in the presence and absence of beam space-charge by a and ctq (both 
measured in degrees per lattice period) [7,8] . The undepressed phase-advance 
(To provides a measure of the strength of the linear applied focusing forces 
of the lattice that is relatively insensitive to the details of the lattice, ctq 
is generally made as large as beam stability will allow - because stronger 
focusing results in smaller beam cross-sectional area leading to smaller, more 
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economical accelerator structures, a can be unambiguously defined by an rms 
equivalent, matched KV equilibrium beam[7,8] where all particles internal 
to the beam have the same phase advance. The ratio ojuQ is a normalized 
measure of relative space-charge strength with ajuQ — > 1 corresponding to 
a warm beam with zero space-charge forces and ojoQ — >■ corresponding 
to a cold beam with maximum space-charge forces. The maximum possible 
current density for a specified beam line-charge density will occur when a ja^ 
is as small as possible. 

Neglecting image charge effects, single particle and beam centroid oscillations 
are stable if ctq < 180°[8]. The parameter space (Tq G (0,180°) and ojoQ e 
(0, 1) can be regarded as potential machine operating points. Envelope models 
predict well understood bands of strong parametric instability when ctq > 90° 
and cr/a"o < 1[7]. The parameter region excluded by envelope instabilities for 
FODO quadrupole transport is indicated (in blue) on Fig. 1. 

Considerations beyond centroid and envelope instabilities exclude further re- 
gions of uq-u parameter space. Transportable current limits based on preser- 
vation of beam statistical emittance and suppression of particle losses for a 
matched beam propagating in a periodic FODO lattice of 84 electric quadrupoles 
were experimentally studied by Tiefenback at LBNL[1,2]. It was found empir- 
ically that transport was stable (i.e., statistical emittance growth and particle 
losses below measurement thresholds) when 

al-a-'<\{\2^^f. (1) 

The additional parameter region this criterion excludes for machine operation 
(partially overlapping the envelope band) is indicated (in red) on Fig. 1. For 
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high space-charge intensity with a /(Jq < 0.5, this hmit is more important than 
the envelope instabihty band because it is encountered first when approaching 
from low ctq- The stability bound (1) has been applied by simply requiring that 
(Jo < 120°/v^ ~ 85°. It is observed that transport becomes more sensitive to 
errors near the boundary of stability. 




(Jo (deg/Period) 



Fig. 1. (color online) Beam stability regions in a FODO quadrupole lattice. 
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2 Pcirticle-in-Cell simulations 

Self-consistent electrostatic particle-in-cell (PIC) simulations have been car- 
ried out for a variety of initial beam distributions launched in a FODO quadrupole 
transport channel with 50% quadrupole occupancy (77 = 1/2) and linear, 
piecewise-constant quadrupole forces. The transverse shce module of the WARP 
code [9] is employed to advance an initial transverse distribution with zero axial 
velocity spread. Applied focusing forces are adjusted for specified ctq. Currents 
are adjusted for specified u/cro using fixed rms emittances {e^ = £y = 50 
mm-mrad). Numerical parameters are set for high resolution (> 100 radial 
zones across the beam core on a square mesh and > 100 residence corrected 
axial steps per lattice period) and good statistics (> 400 particles per cell). 
A cylindrical beam pipe is large enough to suppress particle losses and image 
charge effects. Simulation results are in qualitative agreement with Eq. (1) for 
a wide variety of initial distribution functions. Initial distributions employed 
are rms matched transforms of continuous focusing equilibrium Waterbag and 
thermal distributions [8], KV, or semi-Gaussian distributions. This contrasts 
earlier work where Waterbag and Gaussian loads did not include space-charge 
screening effects and were far from initial force-balance[4] . Only the initial 
KV load employed is a true equilibrium of the periodic focusing channel. No 
exact, smooth equilibrium distributions are presently known for periodic fo- 
cusing channels. 

Parameters to the right of the stability bound [Eq. (1)] and to the left of the 
envelope instability band lead to statistical (rms) emittance growth and parti- 
cle losses. X- and plane average emittance [{e^ + ey)/2] growth can be rapid 
and large as illustrated in Fig. 2(a) for a focusing channel with ctq = 100° 
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and three initial distributions: semi-Gaussian, waterbag "equilibrium", and 
thermal "equilibrium." Much of this emittance growth can be traced to parti- 
cles that evolve significantly outside the beam core as evident from Fig. 2(b) 
which shows the fraction of beam particles which evolve (at one or more 
points) more than 1.25 and 1.5 times the statistical beam edge radius (i.e., 
■sjx'^/r'l + y^jr"^ > 1-25, 1.5 with — 2(a;^)^/^ and — 2(a;^)^/^ calculated 
from the evolving distribution). The similarity of the results for the three 
very different non-equilibrium distributions shows that processes degrading 
the beam are relatively insensitive to the form of the initial distribution in 
deeply unstable parameter regions. Laboratory beams are born off a source 
(injector) and subsequently manipulated to match into a transport channel 
and are unlikely to be any equilibrium form. More detailed analysis of the 
simulation results show that initial beam distortions leading to the statistical 
emittance growth are primarily near the edge of the beam and subsequently 
act to strongly perturb the core. Core perturbations are observed in both the 
local density and temperature profiles. These perturbations typically lack el- 
liptical symmetry and rapidly oscillate into the core with excursions larger 
near the beam edge. Beam envelope matches are not significantly degraded in 
the initial stages of instability. 

A large number of PIC simulations were carried out to better quantify para- 
metric regions of instability. Plane averaged emittance growth contours in ctq 
and a/ao are shown in Fig. 3 for an initial semi-Gaussian distribution. Irreg- 
ular grid simulation points are indicated with dots. All simulations are ad- 
vanced for six undepressed betatron periods, which is sufficient for saturation 
in strongly unstable regimes. Near the stability boundary, emittance growth 
slows and growth factors increase with longer propagation distance. Colors 
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Fig. 2. (color online) PIC simulations of (a) the plane-averaged emittance growth 
for different initial distributions in a FODO quadrupole channel, and (b) the fraction 
of the beam distribution evolving outside the core, (do = 100°, r] = 0.5, Lp = 0.5 
m, fj/fJo = 0.2, and e = 50 mm-mrad). 

show logarithmic scale emittance growth and 1% and 10% threshold contours 
(dashed) are labeled separately. The extent of the envelope instability band 
and Tiefenback's stability threshold are indicated. Results for initial waterbag 
and thermal distributions are similar, but the transition to instability has more 
structure for the waterbag distribution. Strong growth regions in all cases are 
qualitatively consistent with Tiefenback's threshold. Emittance growth can- 
not be attributed to KV-like modes internal to the beam[10,ll]. Much of the 
emittance growth is associated with particles that evolve significantly outside 
the core (see Fig. 2b) rendering any linear internal mode interpretation ques- 
tionable. Also, many KV modes are strongly unstable (instabilities exist for 
a/ao < 0.3985 even in the continuous focusing limit) where no rms emittance 
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growth is observed. KV modes generally predict wrong parametric variations 
of instability (thresholds bend the wrong way). Large internal modes also 
possess little free energy to drive statistical emittance growth[8] and therefore 
may not be dangerous if they saturate at small amplitudes. 
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Fig. 3. (color online) Contours of emittance growth for an initial semi-Gaussian 
distribution in a FODO quadrupole channel. (77 = 0.5, Lp = 0.5 m, and e = 50 
mm-mrad) . 



3 Core-Particle model 




Consider an unbunched beam of ions of charge q and mass m propagating 
with axial velocity jSbC (c is the speed of light in vacuuo) and relativistic 



factor 7b = l/yl — /S^. A linear applied focusing lattice is assumed, self-field 
interactions are electrostatic. Then the transverse orbit x{s) of a beam particle 
evolves according to the paraxial equations of motion [7,8] 

Here, s is the axial coordinate of a beam slice, primes denote derivatives with 
respect to s, and Hix{s) is the linear applied focusing function of the lattice 
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(specific forms can be found in Ref. [7]), and the electrostatic potential is 
related to the number density of beam particles n by the Poisson equation 
= —qn/eo in free-space. Cq is the permittivity of free space. 

The core of the beam is centered on-axis {x — — y), and is uniform density 
within an elliptical cross-section with edge radii rj (henceforth, j ranges over 
both X and y) that obey the KV envelope equations 



' X ' ' y ' j 



Here, Q — gA/(27reom7^/?^c^) — const is the dimensionless perveance, (A = 
qn{x — 0,y — 0)rxry — const is the beam line-charge density), and Sj is the 
rms edge emittance along the j-plane. We take ej = s = const. For a periodic 
focusing channel with lattice period Lp, Kj{s + Lp) — Kj{s), the envelope is 
called matched when it has the periodicity of the lattice, i.e., rj{s + Lp) — 
rj{s). Undepressed particle phase advances are used to set the lattice focusing 
functions Kj using cosctq = (l/2)Tr M where M is the x or y plane transfer 
map of a single particle {Q = 0) through one lattice period. We take the Kj to 
be piecewise constant with occupancy r] e (0, 1] . The matched beam envelope 
flutter varies only weakly with r) but increases strongly with increasing ctq. 
The depressed particle phase advance is calculated as a = e Jo'^ds/rj. 

It can be shown that the flutter of the matched beam envelope for periodic 
FODO quadrupole focusing systems with piecewise constant Kj{s) is given 
approximately (for a/a^ <^ 1) by[12] 

-- - 1 ^ (1 -cosao) 23/2(1 -2,/3)V2 
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Here, rj e (0, 1] is the occupancy of the quadrupoles in the lattice and fx — 
(l/Lp) Jo'^ds Tx- Equation (4) shows that envelope flutter in a quadrupole 
channel depends strongly on ctq and weakly on r] (the variation in rx\xaay./fx in 
77 is ~ 13%). 

For a particle evolving both inside and outside the elliptical beam envelope, 
Eq. (2) can be expressed as 

X -\- KxX — , . X, yo) 



with an analogous equation for the y-plane. Here, Fj are form factors satisfying 
Fj — 1 inside the beam (x^/r^ + y^/r^ < 1) and F^ = {r^ + ry)^Re[S] and 
Fy — —{rx + ry)Y^'m[S] outside the beam {x^/rl -\-y^ /Vy > 1). ^2 is a complex 
variable deflned as S_= ^^21:^2) [1 — ^ 1 — ^-^^^] , where z — x + iy and i — y/^. 

The particle equations of motion (5) are integrated numerically from initial 
conditions. We typically launch particles with initial x and y coordinates out- 
side the beam edge (i.e., x'^/r^ + y'^/Vy > 1) and with initial angles x' and 
y' consistent with coherent flutter motion of core particles extrapolated to 
the location of the particle, i.e., with x' — r'^x/rx and y' — r'yy/ry. Di- 
agnostics include particle trajectories, single particle emittances deflned by 
Ex — \J {x/txY + {xrij. — x'TxY /si {^x = 1 at the core distribution edge), 
stroboscopic Poincare phase space plots, and particle oscillation wavelengths 
calculated from Fourier transforms of orbits. Particle trajectories and phase- 
spaces analyzed in scaled units (e.g., with x-x' projections scaled as xjvx- 
{x'tx — r'^x)/ex) to better illustrate oscillation extents relative to the matched 
beam core. 
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4 Core-Particle simulations 



To illustrate the halo structure, we launch particles along the x-axis of the el- 
liptical beam in specified regions outside the beam edge (e.g., x G [1.1, 1.2]r2.) 
with zero incoherent angle spreads (e.g., x' — r'^x/r^)- Fig. 4 illustrates x-x' 
Poincare phase-spaces for particles launched with x e [1.1, 1.2]rx for fixed ctq 
and two values of cr/cxo: (a) a high value (weak space-charge) well within the 
stable region of Fig. 1, and (b) a low value (strong space charge) in the unsta- 
ble region. The Poincare strobe is one lattice period. Scaled coordinates x/vx 
and {x'r^ — xr'^jzx are plotted to remove envelope flutter. The extent of the 
core is plotted in red. Extrapolations of the range of initial launch conditions 
are indicated in red based on the annular elliptical region formed if the initial 
particle conditions evolved with constant single-particle emittance e^. Note 
the large change in scale between the stable and unstable plots. For the stable 
case, particles diving in and out of the matched envelope remain close to the 
initial launch range and indicate a weak, high-order resonance. For the un- 
stable case, numerous resonances near the core become stronger and overlap 
causing the region immediately outside the core to break up into a stochas- 
tic sea that closely approaches the core. A large, 4-lobe bounding resonance 
(KAM surface) persists that ultimately limits the achievable particle oscilla- 
tion amplitude. The phase advance of particles moving outside the envelope 
is strongly amplitude dependent ranging from a for amplitudes at the core 
boundary to (Tq for very large amplitudes. Strong space charge {cyjaQ <ti 1) 
and large matched envelope oscillations (large ctq) provide a strong pump at 
the lattice frequency. Numerous harmonics of particle orbits near the core 
resonate with the lattice resulting in overlapping resonances that produce a 
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strongly chaotic region that approaches the core. This chaotic sea allows par- 
ticles near the core to rapidly evolve to large amplitudes. 




Fig. 4. (color online) Core-particle Poincare phase-spaces for ctq = 100°, cr/ao = 0.67 
(a), and (t/(To = 0.2 (b). (Lp = 0.5 m, 77 = 0.5, e = 50 mm-mrad). 

A new stability criterion is adopted to estimate where chaotic halo orbits near 
the beam core can degrade transport. When varying cxo and cr/ao, we define the 
stability boundary to be the first point when approached from stable regions 
(low ctq) where particle groups launched near the core (e.g., x G [1.05, l.lOjr^.) 
experience large increases in oscillation amplitude (e.g., Max[x/r2;] increased to 
1.5). Boundary points obtained when particles launched with x/r^ G [1.05, 1.10] 
increase in amplitude by factors of 1.5 (triangles) and 1.4 (squares) are plot- 
ted in Fig. 1. The boundary roughly tracks the region of strong emittance 
growth observed in experiment and simulations until the envelope instability 
band is approached. Results are relatively insensitive to the choice in initial 
group radius and amplitude increase factor. Earlier work by Langiel[13] em- 
ployed a core-particle model to analyze transport limits but implied overly 
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pessimistic stability criteria (ctq < 60° and (yjcjQ > 0.4) seemingly based on 
rough resonance overlap estimates. 




30 60 90 120 150 180 

ctq (degrees/period) 



Fig. 5. (color online) Beam stability boundary calculated from a core-particle model 
for a FODO quadrupole channel. {Lp = 0.5 m, r/ = 0.5, e = 50 mm-mrad). 

Halo properties analyzed persist when particles have finite angular momen- 
tum (not launched on-axis). Particles that leave the core in self-consistent 
PIC simulations generate similar Poincare plots to the core-particle model for 
a variety of initial distributions. Single particle emittance growths of ~ 50 
are possible for particles near the beam edge that enter the halo in unstable 
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regions. If a significant number of near-edge particles enter the halo, this can 
result in strong increases in rms beam emittance and distortions in the beam 
phase-space (both total and core) . Particles leaving the core in unstable regions 
rapidly grow in amplitude over a relatively small number of lattice periods - 
consistent with PIC simulations. Moreover, as observed in simulations and ex- 
periment, this halo induced mechanism for transport degradation is consistent 
with increasing sensitivity to the beam distribution and edge perturbations as 
the threshold region is approached. The core-particle model assumption of a 
uniform density elliptical beam core is reasonable for strong space-charge due 
to Debye screening and phase-mixing of initial perturbations. No periodic, 
nonuniform density equilibria are known and core perturbations are observed 
in PIC simulations to collectively evolve and disperse leaving smaller residual 
fluctuations and a rounded beam edge. Hence the uniform core model can pro- 
vide a good approximation to the average impulse a halo particle experiences 
while traveling through the oscillating core. If the edge of the beam distri- 
bution is not sharp, as is expected for finite (j/o-q, a significant population of 
edge particles can enter the halo and be elevated to large amplitudes in unsta- 
ble regions. Due to envelope flutter, the spatial average temperature of a 
beam with constant emittance Sx will vary as oc sl/rl- Thus, in a matched 
beam envelope the temperature will oscillate with the period of the lattice 
(360° phase advance), decreasing where the envelope contracts and conversely 
decreasing where the envelope expands. For high ctq these fluctuations will 
tend to increase leaving the edge out of force-balance. On the other hand, the 
plasma response of the beam will have characteristic collective phase advance 
(Tp = ^^7^-\/2Q This frequency will generally be incommensurate with and 
much slower than the temperature oscillations showing showing that the edge 
of the beam will have a more difficult time readapting to the focusing kicks 
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as (To increase and envelope flutter becomes larger. 
5 Conclusions 

A core-particle model was used to analyze the previously unexplained ori- 
gin of space-charge related transport limits of beams propagating in periodic 
quadrupole focusing channels. It was shown that when matched beam enve- 
lope oscillations and space-charge strength are both sufficiently large, near- 
edge particles oscillating both inside and outside the matched beam envelope 
become chaotic and can experience large increases in oscillation amplitude. 
This resonance halo need not be tenuous and is distinct from envelope mis- 
match driven halo [14] because the driving oscillation is the fast flutter of the 
matched beam envelope rather than envelope mismatch modes. The matched 
envelope flutter becomes larger with increasing ctq, providing a strong pump 
that further increases as beam space-charge forces become larger. Envelope os- 
cillations also drive large temperature oscillations in the core of the matched 
envelope. Because the collective response of the beam to local force imbalances 
scales with the plasma frequency which is low relative to the lattice frequency, 
it is unhkely that the beam edge can consistently adapt. Lack of edge self- 
consistency in periodically focused beam distributions makes it plausible that 
many near-edge particles can move sufficiently outside the beam core to par- 
take in the resonance. Consequently, large distortions in the beam phase-space 
and large rms emittance growth can result. Stability thresholds based on this 
resonance picture are in rough agreement with experimental measurements 
and simulations. Analogous transport limits to the ones studied here will oc- 
cur in other periodic focusing channels. Generally, unstable parameters will 
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differ due to different scaling of matcfied beam envelope flutter. Envelope mis- 
match also increases driving envelope excursions and introduces additional 
frequencies - likely reducing the region of stable transport. Work is ongoing 
to further clarify the processes described. Further details of this work can be 
obtained on the arXiv e-print server[15] and in future publications. 
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